home *** CD-ROM | disk | FTP | other *** search
/ The X-Philes (2nd Revision) / The X-Philes Number 1 (1995).iso / xphiles / hp48_1 / analemma < prev    next >
Internet Message Format  |  1995-03-31  |  21KB

  1. From: Craig A. Finseth <fin@unet.umn.edu>
  2. Subject:  v04i001:  analemma_cf - Analemma v1.0, Part01/01
  3. Newsgroups: comp.sources.hp48
  4. Organization: University of Minnesota, Networking Services.
  5. Followup-To: comp.sys.hp48
  6. Approved: spell@seq.uncwil.edu
  7.  
  8. Checksum: 1531207743 (verify with brik -cv)
  9. Submitted-by: Craig A. Finseth <fin@unet.umn.edu>
  10. Posting-number: Volume 4, Issue 1
  11. Archive-name: analemma_cf/part01
  12.  
  13. BEGIN_DOC analemma.doc
  14. Craig A. Finseth            fin@unet.umn.edu
  15. 1343 Lafond                Craig.Finseth@nic.mr.net
  16. St Paul MN 55104-2437            +1 612 644 4027
  17. USA
  18. When: 14 January 1992
  19.  
  20. What it does:
  21.  
  22. Many globes show a large "figure 8" somewhere over the Pacific Ocean.
  23. (Oddly, this figure does not show on photographs taken from space...) 
  24. This figure is called an "analemma," and it shows the position of the
  25. Sun (in degrees of declination and minutes of right ascention (that
  26. is, minutes as in "hours and minutes" not "degrees and minutes"))
  27. relative to "where it is supposed to be."
  28.  
  29. For example, on May 14, the Sun is about 3 minutes and 41 seconds
  30. ahead of nominal and about 20 degrees north of the ecliptic.
  31.  
  32. I said "many globes" because some manufacturers are removing this
  33. useful piece of information from their globes. This program makes up
  34. for this loss (as well as providing an excuse to write a program).
  35.  
  36. This program uses the graphics capabilities of the 48S/48SX to draw an
  37. analemma for the year 2000.
  38.  
  39. The equations used are from:
  40.  
  41.     Meeus, Jean (1991) Astronomical Algorithms.  Richmond,
  42.     Virginia: Willmann-Bell.  ISBN 0-943396-35-2.
  43.  
  44. In many cases, the definitions for the functions in this program give
  45. equation numbers from this text.
  46.  
  47. Note: the shape on an analemma varies over time.  Also, an analemma is
  48. *not* symmetric about the origin along either axis.
  49. ------------------------------------------------------------
  50. Basic Operation:
  51.  
  52. Execute "GO".  Wait and watch.
  53.  
  54. When complete, you can enter the graphics environment and use the
  55. cursor and COORD to identify the position in decimal minutes of right
  56. ascention and degrees of declination.
  57.  
  58. The program operates in four phases:
  59.  
  60. 1) Clear the screen and draw the grid (declination lines of -20, -10,
  61. 0, 10, and 20 degrees and a right ascention line of 0 degrees).
  62.  
  63. 2) Draw the analemma curve.  Points are calculated every five days and
  64. a line is drawn connecting the points.
  65.  
  66. 3) Draw the tick marks.  These are drawn the first of every month and
  67. 21 Mar, 20 Jun, 21 Sep, and 20 Dec.  February is assumed to have 28
  68. days.
  69.  
  70. 4) Draw the labels.
  71.  
  72. The program objects all communicate via global variables.  Each
  73. program uses those globals that it needs and stores its result in its
  74. own global variable.
  75.  
  76. Because of these dependencies, the programs must be executed in the
  77. following order:
  78.  
  79.     FIGJD
  80.     FIGT
  81.     FIG\GD
  82.     FIG\Ge
  83.     FIGL0    anytime after FIGT
  84.     FIG\Ga
  85.     FIGE
  86.     FIG\Gd    anytime after FIGT
  87.  
  88. While FIGL0 may be computed anytime after FIGT, FIG\Ga and FIGE
  89. (indirectly) must be computed after FIGL0.
  90.  
  91. As it turned out, this structure made it very easy to test the
  92. modules.  A test Julian Day could be stored in JD, then the other
  93. program calculated and their results examined.
  94.  
  95. The programs can, of course, be used separately.
  96. ------------------------------------------------------------
  97. Accuracy:
  98.  
  99. The program calculates the analemma for the year 2000.  While the
  100. analemma does vary over time, the variation is not much and the
  101. displayed shape wouldn't change for a couple of centuries on each side
  102. of the target year.
  103.  
  104. Most of the programs assume this range of time.  Higher-order terms
  105. are not included in the programs, but could be added with minimal
  106. effort.
  107.  
  108. The horizontal axis is "magnified" by a factor of four.  If you wish
  109. to change this magnification, you will need to edit:
  110.  
  111.     PPAR    to set new limits
  112.     FRAME    to set horizontal endpoints
  113.     LBL    to position the labels
  114.  
  115. If you don't care about the lines and labels, you can drop the
  116. invocations from GO and simply edit PPAR.  Then replot.
  117.  
  118. Internal computations are all performed in decimal degrees.  Where a
  119. computation may result in a number of degrees significantly in excess
  120. of 360, the value is taken mod 360 before being stored.
  121. ------------------------------------------------------------
  122. Object Summary:    
  123.  
  124. GO    runs the program
  125. BASE    computes the basic analemma
  126. E    global variable: result variable for the equation of time
  127. FIGE    computes E
  128. FIGJD    computes the Julian Day
  129. FIGL0    computes L0
  130.  
  131. FIGT    computes T
  132. FIG\Ga    computes \Ga
  133. FIG\GD    computes \GD\Ge and \GDp
  134. FIG\Gd    computes \Gd
  135. FIG\Ge    computes \Ge
  136. FRAME    displays the frame
  137.  
  138. JD    global variable: current Julian Day
  139. L0    global variable: longitude of the Sun
  140. LBL    displays the month labels    
  141. PLACE    computes the coordinates for a day
  142. PPAR    global variable: plotting parameters
  143. T    global variable: current Julian Century
  144.  
  145. TICK    places month tick marks
  146. \Ga    global variable: apparent Solar right ascention
  147. \GD\Ge    global variable: nutation in obliquity
  148. \GDp    global variable: nutation in longitude (p should be psi)
  149. \Gd    global variable: Solar declination
  150. \Ge    global variable: obliquity of the ecliptic
  151. ------------------------------------------------------------
  152. Detailed Interfaces:
  153.  
  154. GO    Stack Input:    none
  155.     Stack Output:    none
  156.     Global Input:    none
  157.     Global Output:    none
  158.     Calls:        BASE, FRAME, LBL, TICK
  159.  
  160.     Saves the flags, sets degrees mode, clears and displays the
  161.     PICT, draws the figure, restores the flags.
  162.  
  163. BASE    Stack Input:    none
  164.     Stack Output:    none
  165.     Global Input:    none
  166.     Global Output:    none
  167.     Calls:        PLACE
  168.  
  169.     Plots the basic analemma.
  170.  
  171. FIGE    Stack Input:    none
  172.     Stack Output:    none
  173.     Global Input:    L0, \Ga, \GDp, \Ge
  174.     Global Output:    E
  175.     Calls:        none
  176.  
  177.     Computes the equation of time (equation 27.1).
  178.  
  179. FIGJD    Stack Input:    day number within year (21 Mar = 0)
  180.     Stack Output:    none
  181.     Global Input:    none
  182.     Global Output:    JD
  183.     Calls:        none
  184.  
  185.     Computes the Julian Day
  186.  
  187. FIGL0    Stack Input:    none
  188.     Stack Output:    none
  189.     Global Input:    T
  190.     Global Output:    L0
  191.     Calls:        none
  192.  
  193.     Computes the Sun's mean longitude (equation 27.2).
  194.  
  195. FIGT    Stack Input:    none
  196.     Stack Output:    none
  197.     Global Input:    JD
  198.     Global Output:    T
  199.     Calls:        none
  200.  
  201.     Computes the Julian Century (equation 21.1).  The origin is
  202.     1 Jan 2000.
  203.  
  204. FIG\Ga    Stack Input:    none
  205.     Stack Output:    none
  206.     Global Input:    L0, T, \Ge
  207.     Global Output:    \Ga
  208.     Calls:        none
  209.  
  210.     Computes the *apparent* Solar right ascention (equations 24.3,
  211.     both unnumbered between 24.4 and 24.5 ("OD" refers to the
  212.     circle with dot varaible), first two unnumbered between 24.5
  213.     and 24.6, 24.6, and 24.8).
  214.  
  215. FIG\GD    Stack Input:    none
  216.     Stack Output:    none
  217.     Global Input:    T
  218.     Global Output:    \GD\Ge, \GDp
  219.     Calls:        none
  220.  
  221.     Computes the nutation in longitude (\GDp, p should be psi) and
  222.     nutation in obliquity (\Ge) (the accuracy of 0.5 minutes
  223.     equations on page 132).
  224.  
  225. FIG\Gd    Stack Input:    none
  226.     Stack Output:    none
  227.     Global Input:    T, \Ge
  228.     Global Output:    \Gd
  229.     Calls:        none
  230.  
  231.     Computes the Solar declination.  Uses simple first order
  232.     approximation.
  233.  
  234. FIG\Ge    Stack Input:    none
  235.     Stack Output:    none
  236.     Global Input:    T
  237.     Global Output:    \Ge
  238.     Calls:        none
  239.  
  240.     Computes the mean obliquity of the ecliptic (equation 21.1).
  241.  
  242. FRAME    Stack Input:    none
  243.     Stack Output:    none
  244.     Global Input:    none
  245.     Global Output:    none
  246.     Calls:        none
  247.  
  248.     Plots the lines of declination (-20, -10, 0, 10, and 20
  249.     degrees) and right ascention (0).
  250.  
  251. LBL    Stack Input:    none
  252.     Stack Output:    none
  253.     Global Input:    none
  254.     Global Output:    none
  255.     Calls:        none
  256.  
  257.     Plots the month name labels.
  258.  
  259. PLACE    Stack Input:    day number within year (21 Mar = 0)
  260.     Stack Output:    complex coordinate of analemma point for that day
  261.     Global Input:    E, \Gd
  262.     Global Output:    none
  263.     Calls:        FIGE, FIGJD, FIGL0, FIGT, FIG\Ga, FIG\GD
  264.             FIG\Gd, FIG\Ge
  265.     computes the coordinates for a day
  266.  
  267. TICK    Stack Input:    none
  268.     Stack Output:    none
  269.     Global Input:    none
  270.     Global Output:    none
  271.     Calls:        PLACE
  272.  
  273.     Plots the tick marks for the first of every month and 21 Mar,
  274.     20 Jun, 21 Sep, and 20 Dec.
  275.  
  276. END_DOC
  277. Checksum: #c6deh
  278. Size: 2859.5
  279. ------------------------------------------------------------
  280.  
  281. BEGIN_RPL analemma.rpl
  282. %%HP: T(3)A(D)F(.);
  283. DIR
  284.   GO
  285.     \<< RCLF DEG
  286. ERASE { # 0h # 0h }
  287. PVIEW FRAME BASE
  288. TICK LBL 7 FREEZE
  289. STOF
  290.     \>>
  291.   BASE
  292.     \<< 0 PLACE 0 365
  293.       FOR d d PLACE
  294. DUP ROT LINE 5
  295.       STEP DROP
  296.     \>>
  297.   E -3.19866112964
  298.   FIGE
  299.     \<< L0 .0057183 -
  300. \Ga - \GDp \Ge COS * +
  301. 'E' STO
  302.     \>>
  303.   FIGJD
  304.     \<< 78 +
  305. 2451545.5 + 'JD'
  306. STO
  307.     \>>
  308.   FIGL0
  309.     \<< T 10 / \-> \Gt
  310.       \<< 280.4664567
  311. 360007.698278 \Gt * +
  312. .03032028 \Gt SQ * +
  313. 360 MOD 'L0' STO
  314.       \>>
  315.     \>>
  316.   FIGT
  317.     \<< JD 2451545 -
  318. 36525 / 'T' STO
  319.     \>>
  320.   FIG\Ga
  321.     \<< 357.5291
  322. 35999.0503 T * +
  323. 360 MOD \-> M
  324.       \<< 1.9146
  325. .004817 T * -
  326. .000014 T SQ * - M
  327. SIN * .019993
  328. .000101 T * - M 2 *
  329. SIN * + .00029 M 3
  330. * SIN * + L0 + \-> OD
  331.         \<< 125.04
  332. 1934.136 T * - 360
  333. MOD \-> \GW
  334.           \<< OD
  335. .00569 - .00478 \GW
  336. SIN * - \-> \Gl
  337.             \<< \Ge
  338. .00256 \GW COS * +
  339. COS \Gl SIN * \Gl COS
  340. SWAP R\->C ARG 360 +
  341. 360 MOD
  342.             \>>
  343.           \>>
  344.         \>>
  345.       \>> '\Ga' STO
  346.     \>>
  347.   FIG\GD
  348.     \<< 125.04452
  349. 1934.136261 T * -
  350. 360 MOD 280.4665
  351. 36000.7698 T * +
  352. 360 MOD 218.3165
  353. 481267.8813 T * +
  354. 360 MOD \-> \GW L LP
  355.       \<< -.00172
  356. HMS\-> \GW SIN *
  357. .000132 HMS\-> 2 L *
  358. SIN * - .000023
  359. HMS\-> 2 LP * SIN * -
  360. .000021 HMS\-> 2 \GW *
  361. SIN * + '\GDp' STO
  362. .00092 HMS\-> \GW COS *
  363. .000057 HMS\-> L 2 *
  364. COS * + .00001 HMS\->
  365. LP 2 * COS * +
  366. .000009 HMS\-> \GW 2 *
  367. COS * - '\GD\Ge' STO
  368.       \>>
  369.     \>>
  370.   FIG\Gd
  371.     \<< T 36525 * 78
  372. - 365 / 360 * 360
  373. MOD SIN \Ge * '\Gd' STO
  374.     \>>
  375.   FIG\Ge
  376.     \<< 23.2621448
  377. HMS\-> .0046815 HMS\->
  378. T * - .000000059
  379. HMS\-> T SQ * - \GD\Ge +
  380. '\Ge' STO
  381.     \>>
  382.   FRAME
  383.     \<< (-51.17,-20)
  384. (51.17,-20) LINE
  385. (-51.17,-10)
  386. (51.17,-10) LINE
  387. (-51.17,0)
  388. (51.17,0) LINE
  389. (-51.17,10)
  390. (51.17,10) LINE
  391. (-51.17,20)
  392. (51.17,20) LINE
  393. (0,-25) (0,25) LINE
  394.     \>>
  395.   JD 2451967.5
  396.   L0 336.902466384
  397.   LBL
  398.     \<< PICT (-32,1)
  399. "Mar" 2 \->GROB GXOR
  400. PICT (-32,16) "Aug"
  401. 2 \->GROB GXOR PICT
  402. (-32,25) "Jul" 1
  403. \->GROB GXOR PICT
  404. (27,25) "Jun" 1
  405. \->GROB GXOR PICT
  406. (27,19) "May" 1
  407. \->GROB GXOR PICT
  408. (27,13) "Apr" 1
  409. \->GROB GXOR PICT
  410. (27,4) "Sep" 2
  411. \->GROB GXOR PICT
  412. (27,-6) "Oct" 2
  413. \->GROB GXOR PICT
  414. (27,-17) "Nov" 1
  415. \->GROB GXOR PICT
  416. (27,-22) "Dec" 1
  417. \->GROB GXOR PICT
  418. (-32,-19) "Jan" 1
  419. \->GROB GXOR PICT
  420. (-32,-9) "Feb" 2
  421. \->GROB GXOR
  422.     \>>
  423.   PLACE
  424.     \<< FIGJD FIGT
  425. FIG\GD FIG\Ge FIGL0
  426. FIG\Ga FIGE FIG\Gd E 4
  427. * \Gd R\->C
  428.     \>>
  429.   PPAR {
  430. (-51.17,-25)
  431. (51.17,25) X 0
  432. (0,0) PARAMETRIC Y
  433. }
  434.   T
  435. 1.15674195756E-2
  436.   TICK
  437.     \<< { 0 10 40 71
  438. 91 101 132 163 184
  439. 193 224 254 274 285
  440. 316 344 } DUP SIZE
  441. \-> d n
  442.       \<< 1 n
  443.         FOR j d j
  444. GET PLACE PIXOFF
  445.         NEXT
  446.       \>>
  447.     \>>
  448.   \Ga 340.091306964
  449.   \GD\Ge
  450. -4.10799632761E-4
  451.   \GDp
  452. -4.47118519728E-3
  453.   \Gd -8.10071364255
  454.   \Ge 23.4387298868
  455. END
  456. END_RPL
  457.  
  458. BEGIN_ASC analemma.asc
  459. %%HP: T(1)A(D)F(.);
  460. "69A20FF75851000000103910339201008688927834320B100010291033920000
  461. 5524631700189B100020B90720339207998279158117449D100020B939203392
  462. 06991672369970149D100010C810339202004696031900430B100040459434B4
  463. 40D9D20E163247A204B2A2339201000000000000010339201000000000000040
  464. 3392010000000000001703392010000000000001903392020000000000010103
  465. 3920200000000000231033920200000000000361033920200000000000481033
  466. 9202000000000003910339202000000000004220339202000000000004520339
  467. 2020000000000047203392020000000000058203392020000000000061303392
  468. 02000000000004430B213078BF18B9C11C432D6E201046D6E2010E6E16329C2A
  469. 2D6E2010E60A132D6E2010A6D6E201046D6E2010A66C7D184E205005C4143454
  470. 443E1C4232EF53293632B21303E100104510339208996575914765110B100040
  471. 050514254047A207792010000000000711591000000000000529779201000000
  472. 000071150100000000000052084E2010854B2A27792000000000000000000000
  473. 0000000000001C6E184E201095B21301A0005005C414345450D9D20E163284E2
  474. 050649474A44484E20406494744584E2040649474B984E20406494743984E205
  475. 0649474C40384E2040649474C884E20406494745484E20406494742984E20105
  476. 4803A2EEDA184E201029E97C193632B2130FB00030C424C430D9D20E1632634E
  477. 17792010000000000002390000000000000010C2A20B0000D41627ED2A2DA5E1
  478. 4E4E1634E17792010000000000002391000000000000610C2A20B0000145776E
  479. D2A2DA5E14E4E1634E17792010000000000002391000000000000520C2A20B00
  480. 00A457C69C2A2DA5E14E4E1634E1779201000000000000720100000000000052
  481. 0C2A20B0000A457E69C2A2DA5E14E4E1634E1779201000000000000720100000
  482. 0000000910C2A20B0000D416979C2A2DA5E14E4E1634E1779201000000000000
  483. 7201000000000000310C2A20B00001407279C2A2DA5E14E4E1634E1779201000
  484. 0000000007200000000000000040C2A20B0000355607ED2A2DA5E14E4E1634E1
  485. 7792010000000000007200000000000000069C2A20B0000F43647ED2A2DA5E14
  486. E4E1634E17792010000000000007201000000000000719C2A20B0000E4F6679C
  487. 2A2DA5E14E4E1634E17792010000000000007201000000000000229C2A20B000
  488. 04456369C2A2DA5E14E4E1634E17792010000000000002391000000000000919
  489. C2A20B0000A416E69C2A2DA5E14E4E1634E17792010000000000002390000000
  490. 000000099C2A20B0000645626ED2A2DA5E14E4E193632B2130A830020C403203
  491. 39202004836642096330D100020A44420339206000000576915420D100050642
  492. 514D45450D9D20E1632779201000000000071159100000000000002977920100
  493. 00000000711501000000000000029893E1779201000000000071159100000000
  494. 00000197792010000000000711501000000000000019893E1779201000000000
  495. 07115900000000000000007792010000000000711500000000000000000893E1
  496. 7792010000000000711591000000000000010779201000000000071150100000
  497. 0000000010893E17792010000000000711591000000000000020779201000000
  498. 0000711501000000000000020893E17792000000000000000001000000000000
  499. 5297792000000000000000001000000000000520893E193632B2130CF1004064
  500. 94743940D9D20E1632339201000008441262320E3FB133920799000000051864
  501. 0E3FB184E201045EEDA190DA1339202990000000000950E3FB184E201045624B
  502. 1EEDA190DA184E2020B93976BA14563284E20103997632DCC0293632B21301C0
  503. 00406494742940D9D20E163284E201045339204000000000525630EEDA133920
  504. 100000000000087090DA133920200000000000563050FA133920200000000000
  505. 0630EEDA1339202000000000000630D4EB1CA4B184E201039EEDA14563284E20
  506. 102997632DCC0293632B21306D00040649474B940D9D20E16323392020000002
  507. 5440521033920300001626314391084E201045EEDA190DA13392020000000000
  508. 00630D4EB133920200000005664082033920400000896700063084E201045EED
  509. A176BA1339202000000000000630D4EB13392020000000561381203392050000
  510. 3188762184084E201045EEDA176BA1339202000000000000630D4EB11C432D6E
  511. 2010D9D6E2010C4D6E2020C405E1632339207990000000002719E3FB1D6E2010
  512. D9CA4B1EEDA1339206990000000002310E3FB1ED2A2D6E2010C4EEDA1CA4B1EE
  513. DA190DA1339205990000000000320E3FB1ED2A2D6E2020C405EEDA1CA4B1EEDA
  514. 190DA1339205990000000000120E3FB1ED2A2D6E2010D9EEDA1CA4B1EEDA176B
  515. A14563284E2020B90797632DCC02339206990000000000290E3FB1D6E2010D95
  516. 05B1EEDA1339205990000000000750E3FB1D6E2010C4ED2A2EEDA1505B1EEDA1
  517. 76BA1339205990000000000010E3FB1D6E2020C405ED2A2EEDA1505B1EEDA176
  518. BA1339204990000000000090E3FB1D6E2010D9ED2A2EEDA1505B1EEDA190DA14
  519. 563284E2020B93997632DCC02EF53293632B2130B430040649474C840D9D20E1
  520. 63233920200000001925753033920400000305099953084E201045EEDA176BA1
  521. 339202000000000000630D4EB11C432D6E2010D4E16323392000000000006419
  522. 1033920799000000007184084E201045EEDA190DA13392059900000000004108
  523. 4E201045624B1EEDA190DA1D6E2010D4CA4B1EEDA13392089900000003999103
  524. 3920699000000000101084E201045EEDA190DA1D6E2010D4ED2A2EEDA1CA4B1E
  525. EDA176BA1339206990000000000920D6E2010D43F2A2EEDA1CA4B1EEDA176BA1
  526. 84E2020C40376BA11C432D6E2020F444E1632339202000000000405210339203
  527. 00000006314391084E201045EEDA190DA1339202000000000000630D4EB11C43
  528. 2D6E2010D9E1632D6E2020F44433920799000000000965090DA1339207990000
  529. 000008740D6E2010D9CA4B1EEDA190DA11C432D6E201069E163284E201039339
  530. 207990000000006520D6E2010D9505B1EEDA176BA1505B1D6E201069CA4B1EED
  531. A1D6E201069505B1DBBF1E97C1BD2B133920200000000000063076BA13392020
  532. 00000000000630D4EB1EF532EF532EF532EF5324563284E2010C897632DCC029
  533. 3632B2130C5300406494744540D9D20E163284E2020A44433920600000005451
  534. 542090DA133920400000000052563050FA14563284E20104597632DCC0293632
  535. B21307700050649474C40350D9D20E163284E201045339201000000000000010
  536. 50FA11C432D6E201099E16323392020000765466408203392050087289670006
  537. 30D6E201099EEDA176BA1339208990000082023030D6E201099624B1EEDA176B
  538. A1339202000000000000630D4EB14563284E2020C40397632DCC02EF53293632
  539. B2130BF00050649474A44450D9D20E163233920100000000000087076BA13392
  540. 0600000055451542076BA14563284E2020A44497632DCC0293632B2130070004
  541. 06494745440D9D20E163284E2020C40333920799000000038175090DA184E201
  542. 0C890DA184E2020B90784E201039505B1EEDA176BA14563284E20105497632DC
  543. C0293632B2130E8000105410339200004692116689139B1000402414355440D9
  544. D20E16324B2A284E205005C41434544B2A23392020000000000056300A132D6E
  545. 201046D6E20104684E205005C414345478BF1E0CF1893E1D13A2083328DBF193
  546. 632B2130690002074F420D9D20E1632916C1993C1F52E147A20E4A2051000000
  547. 0000000000000E4A20510000000000000000000B21300F2E184E2050642514D4
  548. 5484E20402414355484E2040459434B484E2030C424C4743A24A5A1F76C19363
  549. 2B2130ED6C"
  550. END_ASC
  551.  
  552. BYTES: #C6DEh 2863.5
  553.  
  554. BEGIN_UU analemma.uue
  555. begin 644 analemma
  556. M2%!(4#0X+466*O!_A14````!DP$S*1``:(@IAT,CL`$``9(!,RD``%5"-G$`G
  557. M@;D!``*;<`(S*7"9*)=1&'%$V0$``IN3`C,I8)EA)V.9!T'9`0`!C`$S*2``D
  558. M9&DPD0`TL`$`!%1)0TL$G2W@82-T*D`K*C,I$````````!`PDP(!````````G
  559. M!#,I$````````'$PDP(!```````0"3,I(```````$!`PDP("```````R`3,I$
  560. M(```````,!8PDP("``````"$`3,I(```````,!DPDP("```````D`C,I(```4
  561. M````0"4PDP("``````!T`C,I(```````4"@PDP("```````6`S,I(```````3
  562. M0#2P$@.'^X&;',$TTN8"`61M+A#@YF$CR:+2Y@(!;J`QTN8"`6IM+A!`UN8"H
  563. M`6K&UX'D`@503$%#143CP20C_C628R,K,3`>``%4`3,I@)E65QET5A&P`0`$_
  564. M4%!!4@1T*G"7`@$`````<!&5`0``````4))W*1```````!=1$````````"6`-
  565. MY`(!6+2B<I<"`````````````````````,'F@>0"`5DK,1`*``503$%#106=F
  566. M+>!A(T@N4&"4=*1$A.0"!$9)1U1(+D!@E'2TB>0"!$9)1Y-(+E!@E'3$!(/DR
  567. M`@1&24>,2"Y`8)1T5(3D`@1&24>22"X04(0P*NZM@>0"`9*>QY%C(RLQ\`L`4
  568. M`TQ"3`.=+>!A(S;D<9<"`0``````(),``````````2PJL```36%RWJ+26A[DV
  569. MY&%#'G<I$````````#(9````````%L"B`@L`$%1WYBTJK>5!3AXVY'&7`@$`/
  570. M`````""3`0``````4`(L*K```$IU;,FBTEH>Y.1A0QYW*1`````````G$```$
  571. M`````"7`H@(+`*!4YY8L*JWE04X>-N1QEP(!``````!P`@$``````)`!+"JPT
  572. M``!-87G)HM):'N3D84,>=RD0````````)Q`````````3P*(""P`0!">7+"JM9
  573. MY4%.'C;D<9<"`0``````<`(`````````!"PJL```4V5PWJ+26A[DY&%#'G<IV
  574. M$````````"<`````````8,FB`@L`\#1&YRTJK>5!3AXVY'&7`@$``````'`"$
  575. M`0``````<)$L*K```$YO=LFBTEH>Y.1A0QYW*1`````````G$````````"+)\
  576. MH@(+`$!4-I8L*JWE04X>-N1QEP(!```````@DP$``````)"1+"JP``!*86[)+
  577. MHM):'N3D84,>=RD0````````,@D```````"0R:(""P!@5";F+2JMY4%.'CDVX
  578. MLA(#B@,@P`0C,),"`D`X9B20-@,=`""@1"0PDP(&``!09QE%`AT`4&`D%=143
  579. M5-#9`AXV<I<"`0````!P$94!````````DG<I$```````%U$0````````((DYR
  580. M'G<I$```````%U$9````````$'F7`@$`````<!$%`0```````)&8XW&7`@$`'
  581. M````<!&5``````````!W*1```````!=1``````````"`.1YW*1```````!=15
  582. M&0```````!!PEP(!`````'`1!0$````````!F.-QEP(!`````'`1E0$`````U
  583. M```"=RD0```````741`````````@@#D>=RD``````````!`````````E>9<"T
  584. M```````````!``````!0`ICCD6,C*S'`'P`$1DE'DP2=+>!A(S,I$```@$0AA
  585. M)B/@\QLS*7"9````4(%&X/,;2"X00.7>&@FM,9,"D@D`````D`4^OX'D`@%4[
  586. M)K3AWAH)K8'D`@*;DV>K064C2"X0,)EG(\T,DF,C*S$0#``$1DE'D@2=+>!A5
  587. M(T@N$$`UDP($`````"5E`^ZM,9,"`0``````@`<)K3&3`@(``````&4#!:\QT
  588. MDP("``````!@`^ZM,9,"`@``````8`--OL%*&T@N$##IWAI4-H+D`@&2>3;2=
  589. MS"`Y-K(2`]8`0&"4=+1)T-D"'C8RDP("```@100E`3,I,```86(3-!F`Y`(!%
  590. M5.ZMD=`:,RD@````````-M#D&S,I(````%!F!"@PDP($``"8=@!@`T@N$$#E'
  591. MWAIGJS&3`@(``````&`#3;XQDP("````93$8`C,I4```$XAG$DB`Y`(!5.ZM2
  592. M<;8:,RD@````````-M#D&\$TTN8"`9UM+A#`U.8"`DQ0'C8RDP*7"0````!R<
  593. MD3Z_T>8"`9VLM.'>&C,I8)D`````(!/@\QO>HM+F`@%,[JW!2AONK9'0&C,II
  594. M4)D``````"/@\QO>HM+F`@),4.ZMP4H;[JV1T!HS*5"9```````AX/,;WJ+2Z
  595. MY@(!G>ZMP4H;[JUQMAI4-H+D`@*;<'DVTLP@,RE@F0``````DN#S&VTN$-!9+
  596. M4!ONK3&3`I4)`````'`%/K_1Y@(!3-ZBXMX:!;7AWAIGJS&3`I4)```````!S
  597. M/K_1Y@("3%#>HN+>&@6UX=X:9ZLQDP*4"0``````"3Z_T>8"`9W>HN+>&@6U)
  598. MX=X:":U!92-(+B"P.9EG(\T,XE\C.3:R$@-+`T!@E'3$2-#9`AXV,I,"`@``Y
  599. M`)%25P,S*4```#!0D)DU@.0"`53NK7&V&C,I(````````#;0Y!O!--+F`@%-K
  600. M'C8RDP(``````$:1`3,I<)D`````%TB`Y`(!5.ZMD=`:,RE0F0``````%(#DO
  601. M`@%4)K3AWAH)K='F`@%-K+3AWAHS*8"9````,)D9,),"E@D``````0%(+A!`3
  602. MY=X:":W1Y@(!3=ZBXMX:K+3AWAIGJS&3`I8)`````)`";2X0T#0O*NZMP4H;?
  603. M[JUQMAI(+B#`!'.V&L$TTN8"`D]$'C8RDP("``````0E`3,I,````&`3-!F`Y
  604. MY`(!5.ZMD=`:,RD@````````-M#D&\$TTN8"`9T>-M+F`@)/1#,I<)D`````@
  605. MD%:0T!HS*7"9`````(!'T.8"`9VLM.'>&@FM$4PC;2X08.EA(T@N$#`YDP*7J
  606. M"0````!6`FTN$-!94!ONK7&V&@6UT>8"`9:LM.'>&FTN$&!94!N]^^%Y'-NR(
  607. M,9,"`@``````8`-GJS&3`@(``````&`#3;[A7R/^->)?(_XU0F4C2"X0P)AGW
  608. M(\T,DF,C*S'`-0`$1DE'5`2=+>!A(T@N(*!$-),"!@```$4510()K3&3`@0`>
  609. M````)64#!:]!92-(+A!`E6<CS0R28R,K,7`'``5&24=,,`6=+>!A(T@N$$`UC
  610. MDP(!`````````06O$4PC;2X0D.EA(S,I(```9T5F!"@PDP(%@">8=@!@`VTN,
  611. M$)#IWAIGJS&3`I@)```H(`,#;2X0D&E"&^ZM<;8:,RD@````````-M#D&U0VF
  612. M@N0"`DPP>3;2S"#^-9)C(RLQL`\`!49)1TI$!9TMX&$C,RD0````````>'"VA
  613. M&C,I8````%54421PMAI4-H+D`@)*1'DVTLP@.3:R$@-P`$!@E'141-#9`AXV<
  614. M@N0"`DPP,REPF0```#`85Y#0&D@N$,"8T!I(+B"P"8?D`@&3!;7AWAIGJT%EJ
  615. M(T@N$%"49R/-#))C(RLQX`@``44!,RD``&0I$6:8,;D!``1"05-%!)TMX&$CS
  616. MM**"Y`(%4$Q!0T6THC*3`@(``````&4#H#'2Y@(!9&TN$$"&Y`(%4$Q!0T6'V
  617. M^^'`'YCCT3$J@#."O1\Y-K(2`Y8`('#T)-#9`AXVDF$<F</Q)1YT*N"D`A4`6
  618. M``````````#@I`(5````````````L!(#\.*!Y`(%1E)!345(+D`@%#15A.0"^
  619. 8!%1)0TM(+C#`),1T-"JDI?%G'#DVLA(#Q
  620. ``
  621. end
  622. END_UU
  623.